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We discuss the lower bound to the lightest Higgs boson Hi in the minimal supersymmetric ex- 
tension of the standard model (MSSM) with explicit CP violation, and the phenomenology of the 
lightest relic neutralino in the same scenario. In particular, adopting the CPX benchmark scenario, 
we find that the combination of experimental constraints coming from LEP, Thallium Electric Dipole 
Moment (EDM) measurements, quorkonium decays, and B s — > /i/i decay favours a region of the 
parameter space where the mass of Hi is in the range 7 GeV < Mh 1 i$ 10 GeV, while 3 < tan/3 < 5. 
Assuming a departure from the usual GUT relation among gaugino masses (|Mi| <C [M2D, we find 
that through resonant annihilation to Hi a neutralino as light as 2.9 GeV can be a viable dark 
matter candidate in this scenario. We call this the CPX light neutralino scenario, and discuss its 
phenomenology showing that indirect Dark Matter searches are compatible with the present exper- 
imental constraints, as long as m x < Mh 1 /2. On the other hand, part of the range m x Mh ± /2 
which is allowed by cosmology is excluded by antiproton fluxes. 

PACS numbers: 12.60.Jv, 95.35.+d, 14.80.Cp 



I. INTRODUCTION 

Supersymmetry (SUSY) is considered one of the most natural extensions of the Standard Model (SM) , providing 
elegant solutions to puzzles as diverse as the SM hierarchy problem, the high-energy unification of the gauge coupling 
constants and the existence of the Dark Matter (DM) in the Universe. In particular, in R-parity conserving SUSY 
scenarios the lightest neutralino turns out to be an ideal thermal DM candidate providing in a natural way the 
correct amount of Cold Dark Matter (CDM) that is needed to drive structure formation, and which is necessary to 
explain the latest data on the energy budget of the Universe from WMAP Q. 

Unfortunately, our ignorance about the details of the SUSY breaking mechanism implies that phenomenological 
analysis on the SUSY DM depend in general on a huge parameter space with more than 100 independent soft masses 
and couplings. This parameter space is usually drastically resized by making use of as many simplifying assumptions 
as possible. For instance, parameters related to flavour-mixing are supposed to be strongly suppressed to match 
the experimental constraints, and are often neglected. On the other hand, in theoretically motivated setups as in 
specific SUSY-breaking scenarios, like in the minimal Supergravity (SUGRA) [3|, the number of free parameters is 
strongly reduced, improving predictability. However, some other assumptions are simply suggested by simplicity, such 
as taking soft-breaking parameters real. 

In particular, many recent analysis have addressed this latter aspect flflflflflflHEa since it has been 



realized that CP violating phases in the soft terms can considerably enrich the phenomenology without violating 
existing constraints. This is also theoretically motivated by the fact that the smallness of neutrino masses implied by 
observation possibly calls for some exotic source of CP violation additional to Yukawa couplings, in order to explain 
Baryogenesis. 

In the CP-conserving SUGRA scenario, the so-called "stau-coannihilation" , "Higgs-funnel" and "focus point" 
benchmark solutions are well known examples of a situation where, thanks to a combination of different experimental 
constraints, quite simple and well defined phenomenological pictures emerge, at the price of a certain amount of tuning 
3] . In this article we point out that a similar situation occurs in an effective Minimal Supersymmetric extension of the 
Standard Model (MSSM) with all soft parameters fixed at the electro-weak scale, when CP violation and departure 
from unification of gaugino masses are considered. 

In particular, we wish to address here the issues of the lower bound for the mass of the lightest Higgs boson 
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Hi and that of the lightest possible mass for the relic neutralino x, when standard assumptions are made for the 
origin and evolution of its relic density. In fact, by combining present experimental constraints, a very simple picture 
(albeit tuned) emerges, where the mass of the lightest Higgs boson Hi is found to be in the range 7 < Mh 1 < 
7.5 GeV, with the ratio of the two vacuum expectation values almost fixed, tan/3 ~ 3. This range can be relaxed 
to: 7 < Mjji ^ 10 GeV and 3 < tan/3 < 5, with quite mild assumptions on the Thallium EDM. In this scenario, 
resonant annihilations of neutralinos with mass m x ~ through Hi exchange in the s channel can drive their 

thermal relic abundance within the limits coming from observation, for values of m x significantly below those allowed 
in CP-conserving scenarios [13| (to our knowledge, the issue of resonant \ annihilation in the context of CP violation 
was first raised at the qualitative level in Q)- 

In the following we will analyse in detail the implications for direct and indirect DM searches of these light neu- 
tralinos, which we will refer to as the CPX light neutralino scenario, concluding that it is indeed viable scenario, with 
prospects of detection in future experiments. 

The plan of the paper is as follows. In Section II the CPX scenario in the MSSM with CP violation is introduced. 
In Section III we discuss various experimental bounds. Section IV is devoted to the discussion of the cosmological 
relic density of CPX light relic neutralinos, and Section V to their phenomenology in DM searches. Our conclusions 
are contained in Section VI. 



II. MSSM WITH EXPLICIT CP VIOLATION: THE CPX SCENARIO 



The tree level Higgs potential of the MSSM is invariant under CP transformations. However, CP can be explicitly 
broken at the loop level. In the presence of sizable CP phases in the relevant soft SUSY breaking terms, a significant 
mixing between the scalar and pseudo-scalar neutral Higgs bosons can be generated 0, 0, @j. As a consequence 
of this CP-violating mixing, the three neutral MSSM Higgs mass eigenstates, labeled in order of increasing mass as 
Mh 1 < Mh 2 < Mh 3 , have no longer definite CP parities, but become mixtures of CP-even and CP-odd states. In 
this scenario, all masses are usually calculated as a function of the charged Higgs boson mass M H ± , instead of the 
pseudoscalar Higgs mass Ma, which is no longer a physical parameter. Much work has been devoted to studying the 
phenomenological features of this radiative Higgs-sector CP violation in the framework of the MSSM 0, S|- 

Due to the large Yukawa couplings, the CP-violating mixing among the neutral Higgs bosons is dominated by the 
contribution of third-generation squarks and is proportional to the combination: 

3 $sm.(Af (J,) 
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with f — t,b. Here /i is the Higgs-mixing parameter in the superpotential and Af denotes the trilinear soft coupling. 
In particular, the amount of CP violation is enhanced when the product of \Ab t t\ and is larger than the difference 
of the sfermion masses squared. At the two-loop level, also the gluino mass parameter M3 becomes relevant through 
threshold corrections to the top- and bottom-quark Yukawa couplings. This contribution depends on the combination 
Arg(M3 fi) and can be important especially when tan (3 = v<ijv\ is large, where vi and vi are the vacuum expectation 
values of the neutral components of the two Higgs doublets that give masses to up-type and down-type quarks, 
respectively. More CP phases become relevant by including subdominant radiative corrections from other sectors 0] . 

In presence of CP violation, the mixing among neutral Higgs bosons is described by a 3x3 real orthogonal matrix 
O, instead of a 2x2 one. The matrix O relates the electro-weak states to the mass eigenstates as: 

(fa,fo,a) T = 0(H u H 2 ,H 3 ) T . (2) 

We note that the elements 0$ x i and O^i are the CP-even components of the i-th Higgs boson, while O a i is the 
corresponding CP-odd component. 
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The Higgs-boson couplings to the SM and SUSY particles could be modified significantly due to the CP violating 
mixing. Among them, one of the most important ones may be the Higgs-boson coupling to a pair of vector bosons, 
9HiWi which is responsible for the production of Higgs bosons at e + e~ colliders: 

C HVV = gM w (w+W-» + -^-Z^zA J2 g Hi vvHi , (3) 
V Zc w J i=1 

where 

9HiW = cgO^iH- spOfci , (4) 

when normalized to the SM value. Here we have used the following abbreviations: sp = sin/3, cp = cos/3, tp = tan/3, 
etc. We note that the two vector bosons W and Z couple only to the CP-even components O l j >1 2 i of the i-th Higgs 
mass eigenstate, and the relevant couplings may be strongly suppressed when the i-th Higgs boson is mostly CP-odd, 

o 2 ai ~i»ol ii: ol 2l . 

The so called CPX scenario has been defined as a showcase benchmark point for studying CP-violating Higgs- mixing 
phenomena Its parameters are all defined at the electro-weak scale, and are chosen in order to enhance the 

combination in Eq.(TT]). In this scenario, SUSY soft parameters are fixed as follows: 

M Qs = M U 3 = M D 3 = M U = M &3 = M SUSY , 

|/z|=4M SUSY , |A tAT |=2M SUSY , |M 3 | = 1 TeV, (5) 

where, with a usual notation, Q, L, U, D and E indicate chiral supermultiplets corresponding to left- and right- 
handed quarks and leptons. In this scenario tan/3, M H ±, and AfsusY are free parameters. As far as CP phases are 
concerned, we adopt, without loss of generality, the convention Arg(/i) = 0, while we assume a common phase for 
all the Af terms, $a = Arg(A f ) = Arg(Ab) — Arg(A T ). As a consequence of this, we end-up with two free physical 
phases: $,4 and $ 3 = Arg(M 3 ). 

In addition to the parameters fixed by the CPX scenario, we need to fix the gaugino masses M\ t 2 for our study. We 
take them as free parameters independently of M 3 since, for them, we chose to relax the usual relations at the electro- 
weak scale: Mi /Mj = gf / with g^j =gauge coupling constants, which originate from the assumption of gaugino-mass 
unification at the GUT scale. The neutralino \ is defined as usual as the lowest-mass linear superposition of B-ino 
B, IU-ino and of the two Higgsino states H®, H%: 



X = aiB + a 2 tU (3) +a 3 H° + a 4 H°. (6) 

In Ref. [3| it was proved that in a CP-conserving effective MSSM with |Aii| << | A^2 1 light neutralinos of a mass 
as low as 7 GeV are allowed. Indeed, for \Mi\ « |M 2 | the LEP constraints do not apply, and the lower bound on the 
neutralino mass is set by the cosmological bound. As shown in [15], these neutralinos turn out to be mainly £>-inos, 
ai ~ 1 and m x ~ \Mi\, with a small Higgsino component given by: 

M ... R M Z 

- — - ~ sin 0w sm p , (7) 

|ai| /i 

where 9w is the Weinberg angle and Mz is the Z-boson mass. In the following we will assume vanishing phases for 
Mi and M2, and we will fix for definiteness M2=200 GeV (the phenomenology we are interested in is not sensitive 
to these parameters in a significant way). On the other hand, we will vary Mi, which is directly correlated to the 
lightest neutralino mass m x . 

In this work, we rely on CPsuperH 16| for the computation of mass spectra and couplings in the MSSM Higgs 
sector. 
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III. EXPERIMENTAL CONSTRAINTS ON THE CPX SCENARIO 



A. LEP2 searches 



The most relevant feature of the CPX scenario for our analysis is that the lightest Higgs boson H\ can be very light, 
< 10 GeV, with the other two neutral Higgs bosons significantly heavier, Mh 2:3 ^ 100 GeV, when $^ ~ 90° 
and M H ± ~ 130 GeV for moderate values of 3 < tan/3 < 10. In this case, the lightest Higgs boson is mostly CP 
odd and its production at LEP is highly suppressed since Iff^yyl <§C 1 though it is kinematically accessible. On the 
other hand, the second-lightest Higgs H2 can be produced together with a Z boson since its complementary coupling 
9h 2 vv is sizeable. But its mass is close to the kinematical limit ~ 110 GeV and, moreover, it dominantly decays 
into two Hi's. Depending on Mjj 1 , the lightest Higgs boson decays into two b quarks or two r leptons. This leads to 
a dominant production and decay mode containing 6 jets in the final state, a topology which was covered by LEP2 
with a very low efficiency. The similar situation occurs for H1—H2 pair production. Therefore, in the presence of 
CP-violating mixings, a very light Higgs bosons with Mh x ^ 10 GeV could easily escape detection at LEP2. 

For the CPX scenario, taking $a = $3 = 90° and Msusy = 0.5 TeV, the combined searches of the four LEP 
collaborations at ,/s = 91 — 209 GeV reported the following two allowed regions fl?) ]: 

Rl : M Hl < 10 GeV for 3 < tan/3 < 10, 

R2 : 30 GeV < M Hl < 50 GeV for 3 < tan/3 < 10. (8) 

These regions will not be fully covered even at the LHC for tan/3 < 7 (Rl) and < 5 (R2) In our analysis we 

will focus on region Rl. 

We observe that in the scenario analysed by the LEP collaborations one has |/x|=2 TeV. For this large value of 
the neutralino is a very pure J3-ino configuration, with a Higgsino contamination 03 ~0.02 (see Eq. [7]). As will 
be shown in Section IIVI this has important consequences for the phenomenology of relic neutralinos, in particular 
suppressing their annihilation cross section, and restricting the possibility of having a relic abundance in the allowed 
range only to the case of resonant annihilation. So, the exploration of different possibilities with lower values of 
could in principle be very relevant for relic neutralinos. However, this would require a re-analysis of LEP data which 
is beyond the scope of this paper. 



B. Electric Dipole Moments 



CP phases in the MSSM are significantly constrained by the EDM measurements. In particular, the EDM of the 
Thallium atom provides currently the most stringent constraint on the MSSM scenario of our interest. The atomic 
EDM of 205 Tl gets its main contributions from two terms 
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20|: 



dn [ecm] = -585 • d e [ecm] - 8.5 x 1CT 19 [eon] • (C s TeV 2 ) 



(9) 



where d e denotes the electron EDM and Cs is the coefficient of the CP-odd electron-nucleon interaction Cc s — 
Cs ei^e NN . The dots denote sub-dominant contributions from 6-dimensional tensor and higher-dimensional op- 
erators. The above quantity is constrained by the experimental 2 — a upper bound on the Thallium EDM, which 
is Q: 



*Tl\ 



< 1.3 x 10" 24 I 



(10) 



The contributions of the first and second generation phases, e.g. f4 e(i , <&A d . etc., to EDMs can be drastically 
reduced either by assuming these phases sufficiently small, or if the first- and second-generation squarks and sleptons 
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FIG. 1: The Thallium EDM d T \ = dm x 10 24 ecm in the CPX scenario with Msusy = 0.5 TeV in the region M Hl ^15 GeV 
and 3 < tan/3 < 10. 27ie different shaded regions correspond to different ranges of \dr\\, as shown: specifically, the narrow 
region consistent with the current thallium EDM constraint, \drn\ < 1, is denoted by black squares. In the blank unshaded region 
we have |dxi| > 100. The region below the thick solid line is excluded by data on T(IS') decay ]2a l. For comparison, the thin 
line shows an estimation of the same boundary obtained using the tree-level coupling taking O a i=l, i.e. Is^^J = tan/3. Also 

shown are the three contour lines of the rescaled B(B S — > /i/i) = (B s — > fj,fi) x 10 7 : B(B S — » /i/T) = 2 (solid), 20 (dotted), and 
200 (dashed). 



are sufficiently heavy. However, even when the contribution of the first and second generation phases to EDM's is 
suppressed, there are still sizeable contributions to EDMs from Higgs-mediated two loop diagrams [22| . Their explicit 
forms for (d e /e) H and Cg may be found in Ref. 23|, expressed in the conventions and notations of CPsuperH 16|. 

In Fig. [TJ we show the rescaled Thallium EDM dxi = gJti x 10 24 in units of ecm in the Mf^-tan/3 plane. Here, 
we consider only the contributions from the Higgs-mediated two-loop diagrams. Different ranges of |cZti | are shown 
explicitly by different shadings. In the blank unshaded region we have obtained \d,Ti\ > 100. We find that a cancellation 
between the contributions from (d e /e) H and Cs occurs when tan/? < 5. This cancellation is responsible for the narrow 
region denoted by black squares with |dxi| < 1, where it is at the level of about about 5 %. 

Finally, we note that the Thallium EDM constraint can be evaded by assuming cancellations between the two-loop 
contributions considered here and other contributions, such as those from first- and second-generation sfermions 
discussed above. In this way the allowed region shown in Fig. [T] can be enlarged. The amount of cancellation can be 
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directly read-off from Fig. [TJ For instance, in the region |c?ti | < 10 it would be less severe than 1 part in 10 (10 %). 

C. Bottomonium decay 



In the region Rl, see Eq. ([8j), the bottomonium decay channel T(1S) — > jHi is kinematically accessible 24|. 
There are two experimental limits on this process, depending on whether the Hi decays to visible particles [2^ or to 
invisible ones [26[ . The second case is allowed when 2m x < Mh 1 , so that Hi can decay to neutralinos which escape 
detection. On the other hand, in the case 2m x > Mh x also the three-body decay (i.e. with a non-monochromatic 7 



scenario is related to its SM counterpart by 24| : 



spectrum) T(IS') — > "/XX nas been constrained |25j]. The branching ratio for the two-body decay calculated in our 



B(T(1S) -> 1 H 1 ) SUSY =B(T(1S)^ 7 Hi) SM x (g^-J 2 , (11) 

where g ft denotes the Higgs coupling to two 6 quarks given by = — 0«i tan/? at the tree level. This implies 

that the experimental upper bounds on this process can be directly converted to a constraint in the plane tan (3-Mh 1 ■ 
The result is shown in Fig. [TJ where the thin (red) line corresponds to the limit obtained by setting (<?^ g 6 ) 2 = tan 2 j3. 
Finite-threshold corrections induced by the gluino and chargino exchanges can modify the coupling g^^, although this 
effect is negligible at low values of tan/?. Moreover, although for our choice of parameters Hi is mostly pseudoscalar, 
O a i can be smaller than 1 up to about 20%. The thick solid line in Fig. [TJ shows the bottomonium constraint when 
the threshold corrections and O a i are fully included. 

From Fig. [TJ one can see that, when the following constraints are combined: (i) the LEP constraint; (ii) Thallium 
EDM; (iii) the limit from bottomonium decay, the allowed parameter space is reduced to: 

7 GeV < M Hl <7.5 GeV and tan/?- 3. (12) 

This region may be enlarged to 

7 GeV < M Hl < 10 GeV and 3 < tan/3 < 5 , (13) 

if we assume 10 %-level cancellation in the Thallium EDM. 

In light of the above discussion and for defmiteness, from now on we will fix Mh 1 =7.5 GeV and tan/?=3 in our 
analysis. Taking into account the CPX parameter choice of Eq.© with &a = ^3 = 90° and Msusy = 0.5 TeV, this 
implies, in particular: M H ± ~ 147 GeV, Mjj, ~ 108 GeV, M Hs ^ 157 GeV. 

D. Other constraints 

As will be discussed in the following sections, if the pseudoscalar Higgs boson mass is in the range (IT2|) . a CPX light 
neutralino with m x < Mjjj/2 can be a viable DM candidate. Due to their very pure B-'mo composition, and to the 
quite low value of tan /?, neutralinos in this mass range evade constraints coming from accelerators. For instance, in 
the CPX light neutralino mass range the present upper bound to the invisible width of the Z-boson implies |a§ — a|| < 
a few percent, a constraint easily evaded in this case. 

As far as Flavour Changing Neutral Currents (FCNC) are concerned, they strongly depend on the assumptions 
about flavour violation in the squark sector. For instance, assuming squarks diagonal in flavour, the SUSY contribution 
to the b — > 57 decay rate is dominated by chargino-stop and H -W loops, which are strongly suppressed in our case 
by the low value of tan j3 and by the fact that there are no light masses to compensate this. 
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The situation is potentially different in the case of the decay B s — > pp, since its dominant SUSY contribution 
scales as tan 6 (3 \p\ 2 /Mf Ii and may have a resonance enhancement when Hi is so light that Mjf 1 ~ Mb,- Neglecting 
the threshold corrections which are not so important in our case, we estimate the branching ratio based on the 
approximated expression [27 1 



B(B S m) 



I r„ \r H ff, 
64tt 



\c\ 2 {o±+o A al ) 



with 



C 



Gfo. 
V2tt 



V tb V t * 



tan 3 /? 
4 sin 2 9 



w 



Ml 



sin 29; 



A/ 3 



(14) 



(15) 



where A/3 = fs{x2) — fi{xi) with Xi = m|/|/i| and fz{x) — irlogx/(l — x) and 9^ the stop mixing angle [49( . 
In Fig. [H we show three contour lines of the rescaled B(B S — > pp) = (B s — ■> /i/i) x 10 7 : ^(iJs — > /i/i) = 2 (solid), 
20 (dotted), and 200 (dashed). For the parameters chosen by combining the results from LEP2 searches, Thallium 
EDM, and Bottomonium decay, Eq. (fl"2"]) , we get : B(B S — > pp)cpx — 6 x 10~ 7 taking /s s = 0.23 GeV. This is three 
times larger than the present 95 % C.L. limit 28]: B(B S — > /i/i) < 2 x 10~ 7 . This can be easily made consistent with 
the present experimental constraint if some mild cancellation takes place. The "GIM operative point" mechanism 
discussed in Ref. 29( may be an example of such cancellation, when the squark mass matrices are flavour diagonal. 
In particular, we find that B(B S — > fifi)cpx is consistent to the experimental upper bound by choosing 0.8 < p < 0.9, 
where p = to^/Mstjsy is the hierarchy factor introduced in Ref. [29j ]. with rriq the soft mass for squarks of the first 
two generations [501 ] . 

As far as the SUSY contribution to the anomalous magnetic dipole moment of the muon 5a^ USY is concerned, 
as is well known uncertainties in the SM make a comparison with the experimental results difficult. In particular, 
by combining the SM hadronic vacuum polarization results obtained from e + e~ and t + t~ data, the SM calcula- 
tion turns out to be compatible with observation, and the following 2-cr allowed interval for <5afj' 7 " sy is found [15] ]: 



-160 < 8a 



SUSY 



10 11 < 680. The corresponding contribution from light neutralinos in the CPX scenario falls com- 



fortably into this range: 8a^ <jp X x 10 — 1-5 (to estimate this we have assumed for the trilinear coupling of the 
smuon the same value of the trilinear couplings of the third family given in Eq. (J5|)). 



IV. THE RELIC DENSITY 

Taking into account the latest data from the cosmic microwave data (CMB) combined with other observations 2] 
the 2-cr interval for the DM density of the Universe (normalized to the critical density) is: 

0.096 < n m h 2 < 0.122, (16) 

where h is the Hubble parameter expressed in units of 100 km s _1 Mpc -1 . In Eq. lfTB"]) the upper bound on £l m h 2 
establishes a strict upper limit for the abundance of any relic particle. In particular, the neutralino relic abundance 
is given by the usual expression: 



where < o- ann v > = Xf{cr ann w)i n t, (c a i m v )\nt = jj^ (fann v) dT/m x being the integral from the present temperature 
To up to the freeze-out temperature Tf of the thermally averaged product of the annihilation cross-section times the 
relative velocity of a pair of neutralinos <7 ann v, Xf is defined as Xf = ^j- and g+(xf) denotes the relativistic degrees 
of freedom of the thermodynamic bath at xf . For the determination of Xf we adopt a standard procedure [30j . 
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In absence of some resonant effect, the natural scale of the annihilation cross section times velocity a ann v of CPX 
light neutralinos is far too small to keep the relic abundance below the upper bound of Eq. (|16p (in particular they 
are very pure B-inos and their mass is below the threshold for annihilation to bottom quarks, which is usually the 
dominant channel of a ann v for light neutralinos ,15|). However, when m x ~ neutralinos annihilate through 

the resonant channel XX Hi standard particles, bringing the relic abundance down to acceptable values. In 
the Boltzmann approximation the thermal average of the resonant cr ann v to the final state / can be obtained in a 
straightforward way from the following relation among interaction rates: 

-f < o ann v > ies , f =< T( XX -» /) >=< T(xx -» H l )B{H 1 -» /) >= n H J? x 1 Hl) B f , (18) 

where brackets indicate thermal average, T x is the zero-temperature Hi annihilation amplitude to neutralinos and 
the thermal average of this quantity is accounted for by the ratio of modified Bessel functions of the first kind Ki 
and K2, Bf is the Hi branching ratio to final state /, rij = gimfK2(xi)/(2Tr Xi) are the equilibrium densities with 
Xi = rrii/T, T the temperature, and g\ the corresponding internal degrees of freedom, g x = 2, gn 1 =l- The factor of 
1/2 in front of Eq.(JT5J) accounts for the identical initial states in the annihilation. From Eq.lfTg]). and summing over 
final states /, one gets: 



< a ann v > res = TO 5 K f(x") V{Hl)Bx{l Bx)@ 1 ~ " 2 ' ' ! |,,! 



with B x — T x /T{Hi) , T(H{) the total decay amplitude of Hi, while 9 is the Heaveside step function [5lj . 

By making use in Eq. (fT9|) of the approximations K\(z) ~ ^2(2) — (7r/(2z)) 1 / 2 exp(— z), valid for 2> 1, the integral 
over temperature can be done analytically, leading to: 



<a ~ v> „ 4,2 x fnHi) B x (l - B x ) 1 5(6 + 1) 

u ann u ^ms — ^ n 



(20) 



where S = M H J(2m x ), (5 X = y/l - 5~ 2 and erf(x) = ^ exp{-t 2 ) dt. 

The result of our calculation is shown in Fig. [21 where the neutralino relic abundance Q x h 2 is shown as a function 
of the mass m x . In this calculation the annihilation cross section has been calculated including the off-resonance 
contribution to the annihilation cross section. In this way we have checked that, for m x < Mh ± /2, < cr ann v > rcs is 
always the dominant contribution in the calculation of Q, x h 2 , and Ea.(|20[) is an excellent approximation of < a ann v >. 
The off-resonance contribution to the annihilation cross section is also responsible for regularizing the relic density 
for m x > Afffj/2, where < a ann v > rcs is vanishing. The asymmetric shape of the curve in Fig. [2] is due to the fact 
that thermal motion allows neutralinos with m x < Mh ± /2 to reach the center-of-mass energy needed to create the 
resonance, while this is not possible for m x > Mh 1 /2. In the same figure, the two horizontal lines indicate the range of 
Eq. (fT()]) . The point shown as a circle and indicated as "case A" will be used in the following section as a representative 
point to calculate some signals. 

In Fig.[2]the neutralino mass range allowed by cosmology is: 3.15 GeV ~ m x ~ 3.83 GeV. Allowing for the variation 
of Mh 1 within the range of Eqs. (|12ll3|) . this range is enlarged to: 

2.93 GeV < m x < 5 GeV. (21) 

In this scenario the neutralino relic abundance can fall in the range of Eq. (|16j) only with some level of tuning at 
the boundaries of the allowed mass range. For intermediate values of m x either the neutralino is a sub-dominant 
component of the DM, or some non-thermal mechanism for its cosmological density needs to be introduced. Of course 
all our considerations are valid if standard assumptions are made about the evolution of the early Universe (e.g. about 
the reheating temperature at the end of inflation, the energy budget driving Hubble expansion, entropy production, 
etc). 
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3.5 4 

m x (GeV) 

FIG. 2: Relic abundance as a function of the neutralino mass m x for the CPX scenario with Mh 1 =7.5 GeV, tan/3=3, 
Msusy=0. 5 TeV and <E>a = $3 = 90°. The two horizontal lines indicate the interval of Eg. H16V 
Q x h 2 — (f2 m /i 2 ) m in, is discussed as "case A" in Section W\ 



The circle, where 



V. DARK MATTER SEARCHES 

Neutralinos are CDM particles, and are supposed to clusterize at the galactic level. This implies that they can 
provide the DM density which is gravitationally measured in our Galaxy. In particular, in this section we will 
assume for the DM density in the neighborhood of the solar system the reference value po= 0.3 GeV/cm 3 . When 
we compare our calculation of VL x h 2 to the interval of Eq. (fT6|) . we interpret the lower bound on the DM relic density 
(f2 m /i 2 ) m i„=0.096 as the average abundance below which the halo density of a specific CDM constituent has to be 
rescaled as compared to the total CDM halo density. So, whenever £l x h 2 < (fl m h 2 ) m in we assume that neutralinos 
provide a local density p x which is only a fraction of pq. For the determination of the rescaling factor £ = p x / po we 
adopt the standard recipe: 



£ = min[l,Cl x h 2 /(Q m h 2 )min\- 



(22) 



Neutralinos in the halo of our Galaxy can be searched for through direct and indirect methods. In particular, CPX 
light neutralinos are quite hard to detect through direct detection. Direct detection consists in the measurement of the 
elastic scattering of neutralinos off the nuclei of a low-background detector. For the mass range of Eq. (l2"T]) , the most 
stringent upper bound on the neutralino-nucleon coherent elastic cross section cr^""^ 0011 ^ is provided by the CRESST-1 
CT scaiar° n ' ~ 10~ 38 cm 2 52]. This value is much above the cross section expected in our scenario, 



which falls in the range (^^^° )cpx 



experiment [32 

10~ 42 cm 2 . This is due to the fact that the neutralino-Higgs coupling which 
dominates this process is suppressed, since tan/3 is small and x' s are verv P ure -B-inos due to the large value of \p\. 
Moreover, CTg™^ 1 ™ 11 '' is dominated by the exchange of scalar Higgs bosons, while Hi is mostly pseudoscalar. Finally, 
contrary to annihilation, no resonant enhancement is present in the elastic cross section, since scattering proceeds 
through i-channel. 

For this reason in this section we concentrate on the indirect detection of CPX light neutralinos. In our scenario the 
neutralino relic density fl x h 2 is driven below the observational limit by the resonant enhancement of the annihilation 
cross section < a ann v >. The same cross section calculated at present times, < <J ann v >o, enters into the calculation 
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3.4 3.5 3.6 
m x (GeV) 

FIG. 3: Rescaled, zero-temperature neutralino annihilation cross section to two photons, as a function ofm x . For the solid line 
the rescaling factor £ is calculated according to Eq. while for the dashed one £ = 1. The solid horizontal line shows the 

corresponding constraint from a search for a 7 line from the GC from EGRET ]3m. The dashed horizontal line is an estimate 
for the sensitivity of GLAST ]3lA for a similar search, when model A — N2 of Ref. fSql is used to extrapolate HESS data to lower 
energies. Circle: see caption of Fig. O 



of the annihilation rate of neutralinos in our galaxy. This could produce observable signals, like 7's, v's or exotic 
components in Cosmic Rays (CR), like antiprotons, positrons, antideuterons. 

Note, however, that one can have < <r ann v >o^C < <r a nn v >■ I n fact, as already shown in Section [IV! the thermal 
motion in the early Universe (x x ~ xj ~ 20) allows neutralino resonant annihilation when m x < Mh 1 /2. However, 
for the same neutralinos the contribution of the resonance to < a ann v >q can be negligible at present times, since 
their temperature in the halo of our Galaxy is of order x Xj o — 10 -6 Xf. This implies that the annihilation cross 
section can be large enough in the early Universe in order to provide the correct relic abundance, but not so large at 
present times as to drive indirect signals beyond observational limits. 

In the following we will discuss expected signals for 7 rays and antiprotons. The results of our analysis are 
summarized in Figs. I3I4T51 and [6l In all figures observables are calculated in the CPX scenario with Mff 1 =7.5 GeV, 
tan/3=3, AfsusY=0.5 TeV and &a = $3 = 90°, and plotted as a function of the neutralino mass m x . The solid 
lines show our results obtained by adopting the rescaling procedure for the local density explained above, while for 
comparison, dashed lines are calculated assuming p x = po- When rescaling is applied, indirect signals are proportional 
to the combination t;p x < <J ann v >o, so they reach their maximum value when Q, x h 2 = (Cl m h ) m in- 

For all our results we have adopted, as a reference model, a Navarro-Frenk- White (NFW) profile for the DM density: 



P{ } ~ Pa (r/a)7[l + (r/a)«]0J-r)/° ' ^ J 

where ro = 8.5 kpc is the distance of the solar system from the Galactic Center (GC), a — 25 kpc is the core radius, 
while (a,f3, 7) = (1,3,1). 
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FIG. 4: Continuum 7 flux from the GC calculated for E-,=122 MeV. For the solid line the neutralino local density is rescaled 
according to Eq. 122\) . while for the dashed one, f = 1. The fluxes are compared to the corresponding measurement from EGRET 
fs^j ], shown as the horizontal solid line. Circle: see caption of Fig. [£| 



A. Gamma flux from the Galactic Center 

In Figs. [3] and [5] we show our results for a 7 signal from the GC, which is the most promising source of 7' from 
neutralino annihilation. Note that signals from the GC are proportional to the line-of-sight integral (i.e., performed 
on a path pointing from the observer to the 7 source) J = 1/AJ7 j l p 2 dl d£l, where f2 is the pointing angle 
of observation in the sky and Af2 is the experimental angular resolution. The quantity J is very sensitive to the 
particular choice of density profile, and may span several orders of magnitude, especially for those models that diverge 



in the origin, as for the NFW, where a cut-off radius r cut is needed [40j. In particular, in our calculation we use: 
r C ut = pc. 

The 7 signal from neutralinos takes two contributions: a line with E 7 = m Xl produced by direct annihilation of 
%'s to two 7's, and a continuum, which is mainly due to the annihilation of 7r°'s produced in the fragmentation and 
decay of other final states (quarks, gluons and r's). 

In the first case, the zero-temperature annihilation cross section of neutralinos to photons < cr ann v >o,77 is usually 
suppressed, since it takes place at the one-loop level. However, in our case the contribution of the Hi resonance, 
XX ~ * H\ — > 77 can lead to a strong enhancement of the signal, as is evident from Fig. [3j where our calculation of 



£ 2 < a ann v >o,77 is compared to the upper bound (horizontal solid line) on the same quantity from EGRET 36 1 



(such analysis has been performed on a region of 10°xl0° around the GC, which, for our assumptions, implies J ~ 
120 GeV 2 cm -6 kpc). In this figure the solid line shows our result when £ 2 is calculated according to Eq. ([22"|) . while 
for the dashed line £ = 1. 

In the same figure we also show with a horizontal dashed line an estimate for the prospect of detection of the 
same quantity with GLAST (in this case we have assumed an angular resolution of A = 1CP J sr in the calculation 



of J, leading to J ~ 2500). This estimate is somewhat uncertain, since HESS 4l| has detected a TeV source of 



7's in the GC which is likely to be of standard origin, representing a background for DM searches, and potentially 



making detection of new physics in that region more difficult 381 ] . Our estimate of the background is obtained by 



extrapolating the HESS source to lower energies by using the model A — N2 described in 38| (we have made the 



same assumptions also to estimate the horizontal solid line in Fig. [6]). Our choice of model A — N2 is optimistic, 
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FIG. 5: Left: x 2 calculated for the top-of-atmosphere p flux compared to the experimental data from BESS and AMS JIH] 
(explicitly shown in the right-hand panel). For the solid line the neutralino local density is rescaled as explained in SectionW\ 
For the dashed one, £ = 1. The solid horizontal line indicates the 99.5 % C.L. upper bound for the \ 2 ■ The comparison of the 
\ 2 with its upper limit allows to set an upper bound to m x , indicated by case B and shown as a triangle. On the other hand 
the circle indicates case A, introduced in Fig. [J[R.ight: Top-of-atmosphere antiproton flux as a function of the kinetic TOA 
energy of p's, for the cases A (solid curve) and B (dashes), shown in the left panel as a circle and a triangle, respectively. In 
both cases the lower curve shows the contribution from primaries produced by neutralino annihilation, while the upper curve is 
the total flux where the primary contribution is added to the standard secondary one. The 32 experimental points JIdl l are the 
same that are used to calculate the x 2 shown in the left-hand panel: full circles: BESS 1995-97; open squares: BESS 1998; 
stars: AMS. 



since it implies the smallest extrapolated background at low energies among those discussed in 38]. Of course, a 
more conservative choice for the model adopted to explain the HESS source could make prospects of DM detection 
for GLAST in the GC much worse. 

As far as the continuum signal is concerned, we have calculated the 7 yield from the final states of neutralino 
annihilation using PYTHIA [^J . The result of the calculation is shown in Fig. 21 where we have evaluated the 7 flux 



6 7 from the GC for E 1 — 122 MeV, and compared it with the corresponding flux measurement from EGRET 39j |. 



shown as a horizontal solid line (this particular energy bin is within the range where the data are well fitted by a 
standard background. In this case the flux has been measured in an angular region of 10° x 4° around the GC, which 
corresponds, for our choice of parameters, to J ~ 184 GeV 2 cm~ 6 kpc). 

From Figs. 13141 we can conclude that, with reasonable choices for the DM density profile, 7 signals in our scenario 
are compatible with observations. In both figures we have indicated with a circle the point indicated as "case A" in 
Fig.H 



B. Antiproton flux 



As far as light ncutralinos are concerned, a particularly stringent limit is provided by the flux of primary antiprotons 
that are produced from the hadronization of neutralino-annihilation final states 43j, [4J, |45( . This is due to the fact 



that the p flux observed experimentally is quite in agreement with that expected from p secondary production from 
so that not much room is left for exotic contributions. Moreover, as all other annihilation processes, 



cosmic rays 



the primary neutralino signal scales with the neutralino number density tx so it is enhanced for light masses. 

However, once they are produced in the DM halo, primary p's interact with the magnetic field of the Galaxy, and a 
complex propagation model is needed in order to calculate the fraction of them that reaches the Earth. Unfortunately, 
the main parameters of the propagation model are fixed by using secondary CR data (such as the B/C ratio) which 
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for the cases A (solid curve) and B (dashes) indicated in Fig. [3 arwi|^| as a circle and a triangle, respectively. The horizontal 
line is an estimate of the background for GLAST (see text). 



mainly depend on the galactic disk, while primary p's from neutralino annihilation are produced in the galactic halo. 
This induces uncertainties in the primary flux as large as two orders of magnitude 44( . In particular, for our analysis 
we have used the public code provided by Ref. 45] for the p propagation. Although the code of Ref. [45} is a simplified 
one, where in particular energy-redistribution effects are neglected, it serves well our needs for checking the viability 
of our scenario. In the left panel of Fig. [5] we show the x 2 calculated by comparing the sum of the primary and 
secondary top-of-atmosphere (TOA) p fluxes to the experimental data from BESS and AMS. In particular, for the 
calculation of the \ 2 we have used the same 32 data points from BESS 1995-97, BESS 1998 and AMS 1998 Q shown 
in the right panel. For our calculation we have used a solar modulation parameter <j>= 500 MV, corresponding to the 
period of minimal solar activity when these experiments have taken data. Moreover, in order to be conservative, we 
have used the minimal propagation model in Table III of Ref. 45j, i.e. 5=0.85, iio=0.0016 kpc 2 Myr -1 , L=l kpc, 
V c =13.5 km/s and Va — 22.4 km/s for the propagation parameters S and Kq, for the size of the diffusion zone L, 
and for the galactic wind V c . In absence of a SUSY contribution we find x 2 —30, which confirms the good agreement 
between the data and the standard secondary production (for this latter flux we use the quantity calculated in [45]]). 
As a conservative upper bound for the \ 2 we take \ 2 = 60, which for 32 degrees of freedom implies a statistical 
disagreement at the level of 99.5 %. This value is shown in Fig. as a solid horizontal line. 

From Fig. [5] one can see that, even when the rescaled neutralino local density is used, the \ 2 from p data exceeds 
60 for m x ~ 3.77 GeV. This values for the neutralino mass is indicated with a triangle (case B), along with the circle 
which indicates case A introduced in Fig. [2] Both cases are shown in the right-hand panel, where the upper curves 
show the total p fluxes, while the lower curves show the SUSY contributions. 

From the discussion of this section, where rriH 1 =7.5 GeV and tan/3=3, we obtain for m x the allowed range: 3.15 
GeV < m x < 3.77 GeV. The boundaries of this range correspond to the cases A and B introduced previously. 
Assuming for Mh 1 and tan/3 the range of Eqs. (TT21 [T3")) this interval for m x is enlarged to 2.93 GeV < m x < 5 GeV. 

We conclude this section by showing and example for the prospects of detection for cases A and B in future DM 
searches: in Fig. |6] the 7 continuum flux from the GC is estimated for GLAST (we have assumed here an angular 
resolution of A = 10~ 5 sr in the calculation of J, leading to J ~ 2500). This flux is compared to an estimate of 
the background for the same detector, shown as a horizontal line, calculated with the same assumptions as in Fig. [31 
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where model A — N2 of Ref. 



is used to extrapolate HESS data to lower energies. 



VI. CONCLUSIONS 



In the present paper we have discussed the lower bound to the lightest Higgs boson Hi in the MSSM with explicit CP 
violation, and the phenomenology of the lightest relic neutralino in the same scenario. In particular, we have examined 
the parameter space region Ms 1 < 10 GeV and 3 < tan/3 < 10, in the CPX scenario, an interval which has not been 
excluded by the combined searches of the four LEP collaborations. We find that the combination of experimental 
constraints coming from Thallium EDM measurements and quorkonium decays restricts the region allowed by LEP 
to: 7 GeV < Mh 1 ^5 7.5 GeV, tan/3 ~ 3. In this range, the branching ratio B(B S — > is compatible to the present 
experimental upper bounds provided some moderate cancellation is allowed between the stop-loop contribution and 
that of other squarks. Furthermore, the allowed parameter space can be relaxed to 7 GeV < Mh 1 ^ 10 GeV and 
3 < tan (3 < 5 if a cancellation less severe than 1 part in 10 is also assumed in Thallium EDM between two-loop 
contributions and, for example, those depending from first- and second-generation phases. 

For the above choice of parameters and assuming a departure from the usual GUT relation among gaugino masses 
(\Mi\ <§; I-M2I) we find that neutralinos with 2.9 GeV < m x < 5 GeV can be viable DM candidates. We refer to 
them as CPX light neutralinos. In particular, in the CPX scenario the neutralino is a very pure J3-ino configuration, 
suppressing its Higgs-mediated cross sections. However, through resonant annihilation to Hi the thermal relic density 
of neutralinos with m x ~ Mf^/2 can be either tuned within the range compatible to WMAP or virtually erased, 
allowing for alternative non-thermal mechanisms. The cosmologically allowed range for m x extends to m\ < 
due to the effect of thermal motion in the early Universe. We also have discussed the phenomenology of CPX light 
neutralinos, showing that signals for indirect Dark Matter searches are compatible with the present experimental 
constraints, as long as m x < Mjj 1 /2. On the other hand, part of the range m x > Mji x j1 allowed by cosmology is 
excluded by antiproton fluxes. 

Finally, we note that our study shares some generic features with other models such as the next-to-minimal 
supersymmetric model (NMSSM), in which a light Higgs boson may escape LEP searches and the observed amount 
of DM is explained through neutralino resonant annihilation 
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